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By means of the molecular dynamics simulation on gradual cooling processes, we inves- 
tigate magnetic properties of classical spin systems only with the magnetic dipole-dipolc 
interaction, which we call dipolar systems. Focusing on their finite-size effect, particularly 
their boundary geometry dependence, we study two finite dipolar squares cut out from a 
square lattice with $ = and 7r/4, where $ is an angle between the direction of the lat- 
tice axis and that of the square boundary. Distinctly different results are obtained in the 
two dipolar squares. In the $ = square, the "from-edge-to-interior freezing" of spins is 
observed. Its ground state has a multi-domain structure whose domains consist of the two 
among infinitely (continuously) degenerated Luttinger-Tisza (LT) ground-state orders on a 
bulk square lattice, i.e., the two antiferromagnetically aligned ferromagnetic chains (af-FMC) 
orders directed in parallel to the two lattice axes. In the 3? = 7r/4 square, on the other hand, 
the freezing starts from the interior of the square, and its ground state is nearly in a single 
domain with one of the two af-FMC orders. These geometry effects are argued to originate 
from the anisotropic nature of the dipole-dipole interaction which depends on the relative 
direction of sites in a real space of the interacting spins. 



1. Introduction 

In recent years, systems consisting of arrayed single-domain ferromagnetic nanoparticles 
have attracted much attention as a possible element with high storage density. 1 They exhibit 
rich variety of magnetic phenomena caused by the interplay between the dipole-dipole interac- 
tions and the magnetic anisotropic energy of magnetic moments of nanoparticles. 2 There also 
appeared the works on finite-size effects on such systems. 3-7 In order to get a deeper insight 
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into these rather new and complicated phenomena, we have been studying magnetic proper- 
ties of the systems of magnetic moments between which only the dipole-dipole interaction is 
present, which we call dipolar systems. Particularly, we have focused on dipolar systems of a 
finite size, and have found a peculiar "from-edge-to-interior freezing" in a dipolar cube by a 
finite-temperature molecular dynamics (MD) simulation. 6 As the temperature is decreased, 
magnetic moments on each edge of the cube predominantly start to freeze and then domains 
with short-raged orders corresponding to the bulk Luttinger-Tisza (LT) orders 8 grow from 
edges to the interior. This phenomenon is a very novel one, which, to our knowledges, has not 
been observed in uniform magnetic systems with ordinary short-ranged exchange interactions. 

The magnetic dipole-dipole interaction is quite ubiquitous and has been well known for 
long. It is long-ranged, i.e., the strength of the interaction between magnetic moments (which 
we call simply as spins hereafter) with relative distance r is proportional to r~ 3 . It has more 
peculiar nature though not so emphasized in general. Namely, not only the strength but 
also the sign of the interaction depends on the direction of relative positions in a real space 
of the interacting spins. This anisotropic nature intrinsically gives rise to the frustration 
effect even if spins lie periodically on a certain lattice, and it yields different ground-state 
magnetic orders depending on the lattice structure as pointed out by LT more than a half 
century ago. 8 [Recently, the lower energy states than the LT ground state have been found 
for the BCC dipolar lattice. 9 ] We can then naturally expect that this anisotropic nature, 
combined with the geometry (boundary) effect, plays a key role on determining the magnetic 
properties of finite dipolar systems. Actually, a variety of phenomena in various finite dipolar 
systems have been found in our preliminary study. 10 There have appeared, on the other hand, 
many numerical works which point out the important role of the dipole-dipole interaction 
on the magnetic properties of an individual ferromagnet where the ferromagnetic exchange 
interaction dominates the dipole-dipole interaction. 2, 11-18 However, the proper understanding 
of genuine dipolar systems of a finite size is quite interesting by itself, and is also considered 
to be one of the fundamental bases for technological development of nanomagnetism. 

In the present paper, we study the ground-state and the freezing characteristics of finite 
dipolar squares cut out from a square lattice, which are the simplest finite dipolar systems. The 
corresponding bulk (infinite) dipolar square lattice has the continuous 0(2) degeneracy in its 
ground state. 19 ' 20 The ground-state order is a 2-dimensional (2D) version of the LT order on 
a simple cubic lattice. 8 It is characterized by the order parameter A4 and the state parameter 
0, where A4 is the magnitude of the order parameter vector M whose a-th component is 
defined by 
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6 = arctanf— -r-). 



(2) 



In Eq. (1), Sf and rf are respectively the a-th component and the coordinate in unit of the 
lattice distance of the 2-th classical Heisenberg spin. The definition of Eq. (2) is introduced 
by taking into account the fact that, in the ground state, all the spins are in parallel with 
the lattice plane which is specified as x-y plain. Three examples of the LT order are shown 
in Figs. 1.1 (a-c). We call the LT orders with O = (a) and ir/2 (b) the antiferromagnetically 
aligned ferromagnetic chains (af-FMC) order, and those with G = 7r/4 (c) and 37r/4 as the 
micro-vortex order. 

A main purpose of the present work is to clarify the role of the unique anisotropic nature 
of the dipole-dipole interaction mentioned above on determining magnetic properties of square 
dipolar systems of a finite size. For this purpose, with focusing on their boundary-geometry 
dependence, we study two types of square cut out from a square lattice with $ = and 7r/4 
with $ being an angle between an edge of the square and the lattice axis (see Fig. 1.1(d) 
and (e)). By a similar MD simulation to the one we performed previously, 6 we have found 
significantly different ground states in the two dipolar squares as schematically shown in 
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Fig. 1.1. Spin configurations in the ground state of the bulk dipolar square lattice (a-c). They have 
the 2x2 periodicity, i.e., the LT order which have the continuous 0(2) degeneracy: = (a), 
7r/2 (b), and ir/4 (c). Dipolar squares cut out from the lattice with $ = and $ = 7r/4 are shown 
in (d) and (e), respectively. 
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Fig. 1.2. Schematic pictures of spin configurations in the ground state for $ = (a), 7r/4 (b). The 
arrows denote the spin direction. 

Fig. 1.2. 

The ground state of the $ = dipolar square (Fig. 1.2(a)) consists of four domains with 
the two af-FMC orders shown in Figs. 1.1(a) and (b). The spins of each af-FMC order align 
in parallel to the corresponding lattice axis. This af-FMC order is attributed to the symmetry 
reduction of the spin Hamiltonian from global 0(2) in the corresponding bulk system to Z2 on 
four edges of the square. The from-edge-to-interior freezing starts nearly from a temperature, 
denoted as T* , where the LT short-range order becomes comparable with the linear dimension 
of the system L. These results are similar to, and naturally expected from our previous results 
on the dipolar cube. 6 The ground state of the $ = 7r/4 dipolar square (Fig. 1.2(b)), on 
the other hand, consists of a single domain with small modulations started from the top 
and bottom corners (not shown in the schematic figure). Each domain has an af-FMC order 
along one of the diagonals of the square, which is one of the lattice axes. This single domain 
configuration comes out from an almost complete disappearance of the boundary effect from 
the zig-zag edges. The symmetry reduction to Z2 occurs only around each corner and it 
stabilizes the ground state. The difference of the symmetry reduction on edges between the 
two squares is reflected in their freezing characteristics which differ significantly from each 
other. 

The present paper is organized as follows. In §2 we explain our model and numerical 
method. After pointing out common thermodynamic behavior of the two dipolar squares in 
§3.1, we discuss properties of the ground state and a few low lying states and the freezing 
characteristics of the $ = and = 7r/4 dipolar squares in §3.2 and §3.3, respectively. Section 
4 is devoted to summary and discussion of the present work. 
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2. Model and Method 

The model we study consists of classical Heisenberg spins, {Si}, which are arrayed in a 
finite square. They correspond to magnetic moments of single-domain ferromagnetic particles 
normalized by magnetization M s of the latter. The dipole-dipole interaction energy between 
spins is written as 

Si- 1 -**®* (3) 

<i,j> »i 

Here, rjj denotes the dimensionless length between sites i and j normalized by the lattice 
constant, a, and e^,- the unit vector along the direction from site i to site j. The value of M 2 



,3 



is contained in the coupling constant, J, and J/a is set as the energy (and temperature with 
fee = 1) unit of the present model. For example, in a Ni nanoparticle array with a 30 nm 
particle diameter and a ~ 100 nm, 21 J/a 3 is estimated to be about 5 K. If Co is replaced for 
Ni, J/a 3 is expected to become of the order of several tens of K. 

The freezing characteristics of the model are analyzed using the Landau-Lifshitz-Gilbert 
(LLG) equation, 22 

Axki + F^U. (4) 



dt 1 + a 2 * 1 ' M 



with the effective magnetic field, ifeff.t) given by 

H eSi = - — ?lL + f u (5) 
M s dSi 1 ^ 

where 7i is the dipole-dipole interaction of Eq. (3). To include the heat-bath effect of temper- 
ature T, we introduce random forces {fi} which obey the distribution laws, 

f) = 0, (fi(t) ® fj(0)) = ^TaMA^t)!, (6) 

where a and 7 are the so-called Gilbert damping constant independent of T and the gyro- 
magnetic constant, respectively. 

In the present study, we solve the above set of equations by using the Euler scheme. The 
typical Larmor period, a 3 M s /jJ, multiplied by (1 + a 2 ), is used as a time unit. We further set 
the damping constant a to 0.18, and the time step of integrating Eq. (4), At, to 0.0088. With 
these values of the parameters, the Larmor precession of a spin due to the internal field of 
averaged magnitude damps in a time comparable to one period of the precession (~ 200At) at 
T = 0. We perform typically 8 cooling runs with a fixed cooling rate, in which we use different 
sets of random numbers generating their initial configurations as well as fluctuation forces. 
For each cooling run, the temperature is initially set to about 1 (= J/a 3 ) and is decreased by 
a step of AT = 0.0015. At each temperature, Eq. (4) is integrated over a period of r = 6 x 10 4 
At. The quantities of interest Q at each temperature are time averaged over the whole period 
of r, and are denoted by (Q). We regard it a quasi-equilibrium average of Q in a finite dipolar 
system, postponing the reasoning of our chosen value of r to §4. The obtained averages are 
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Fig. 3.1. Specific heat C of the dipolar squares as a function of temperature. 

almost independent of the cooling process at high temperatures but not necessarily at low 
temperatures due to the intrinsic frustration effect of the dipole-dipole interactions. We call 
an average of (Q) over different cooling runs as the thermal average and denote it as [(Q)] r - 

In order to investigate the boundary effect, we carry out simulations on $ = dipolar 
squares whose linear size along an edge, L e , are set to 15 and 16, and on $ = 7r/4 dipolar 
squares whose linear size along a diagonal direction, L^, are set to 22 and 23. The total 
number of sites, N, of these squares are comparable: ./V = 225, 256 for $ = respectively 
with L e = 15, 16, and N = 264, 265 for $ = vr/4 respectively with L d = 22, 23. The $ = 
square with L e = 15 is enough large for the from-edge-to-interior freezing to be observed. 

Before going into detailed discussions of the simulated results, we note here that the 
magnitudes of fluctuations, ({Sf — (Sf)) 2 ), are rather isotropic, though spins lie almost within 
a plane of the square on average, i.e., (Sf) ~ at all temperatures we have examined. In 
the present work we do not explicitly examine roles of fluctuation of Sf, or in other words, 
differences between magnetic properties of the Heisenberg and XY dipolar models, leaving 
them future problem. 

3. Results 

3.1 Thermodynamic quantities 

Let us start our discussion with the thermal average of the specific heat, C, obtained by 
the present simulation. It is shown in Fig. 3.1. Although the data of <3? = dipolar squares are 
rather fluctuating, the four sets of data behave similarly and exhibit a peak at temperature 
T = T* ~ 0.50 ~ 0.55. The latter is interpreted as a temperature, at which the averaged 
length scale of the LT short-ranged orders becomes comparable to the system size. We may 
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Fig. 3.2. The magnitudes of freezing parameter S as a function of T. 

say that C as well as T* do not depend sensitively on the boundary geometry. The freezing 
parameter S, on the other hand, behaves differently in the $ = and <!> = tt/4 dipolar squares 
as seen in Fig. 3.2. Here we define S as 

S=^I>]r, (7) 

i 

Si = |$)| with 0i) = \J Si{t)dt. (8) 

Concerned with the order parameter M. similarly defined as 

M = [M LT ] T with M LT = -\(M)\, (9) 

with M being the order parameter vector of Eq. (1), Mlt turns out to depend not only on the 
geometry of the squares but also on a cooling run for a fixed geometry (see Fig. 3.3 below). 
These results imply the difference in the freezing characteristics and so the properties of the 
ground state as well as low-lying excited states of the two dipolar squares whose details we 
discuss separately below. 
3.2 $ = dipolar squares 

We show the order parameter Mlt of these dipolar squares observed in eight independent 
cooling runs in Fig. 3.3. For the L e = 15, (16) square there appear four (two) branches at 
low temperatures. The spin configurations at the lowest temperature for these branches are 
shown in Fig. 3.4. The ground-state configuration of the L e = 16 (even L e ) square shown in 
Fig. 3.4(a-i) corresponds to the schematic one shown in Fig. 1.2(a): spins on four edges and on 
four corners draw a closed loop. An almost identical pattern appears on a surface of the dipolar 
cube examined before. 6 In the excited state shown in Fig. 3.4(a-ii), the direction of the af-FMC 
order in one of the four domains is reversed and there appear two corner spins directed to the 
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Fig. 3.3. The order parameter Mlt of $ = squares in each cooling run: (a) L e = 16 and (b) 
L e = 15. 



diagonals of the square, which we call diagonally-directed (d-d) spins hereafter. In the L e = 15 
(odd L e ) square, in which the af-FMC order having the 2x2 periodicity mismatches with L e , 
we have observed four states at lowest temperatures. The ground state is characterized by 
two d-d spins on one of the diagonals as seen in Fig. 3.4(b-i), the lowest excited state by two 
d-d spins on one of the edges (b-ii), and the further excited state by no d-d spins but with 
a point defect (b-iii) or a dislocation (b-iv) of the af-FMC order inside the square. In both 
even and odd squares, the difference between the total energies of the ground state and the 
low-lying excited state(s) is only of the order of L°. 

In order to understand the energetics of the low-lying energy states simulated, it is easier 
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Fig. 3.4. Spin configurations at the lowest temperature obtained by the cooling processes shown in 
Fig. 3.3: (a) L e = 16 and (b) L e = 15 



for us to consider first a set of states, in which all spins on the square align as a single domain 
of the perfect LT order with a fixed value of 0. We call them the 0-LT states and denote their 
local energy at each site i as ej(0). The sum of the latter over i is twice of the total energy 
(Eq.(3)) of the 0-LT states and is denoted as £(&)■ All {ej(0)}'s are shown in Fig. 3.5. 
Reflecting the 0{2) symmetry of the LT order combined with the square geometry of the 
system, ej(0) has the periodicity written as 

A £i (0) = £i (0) - €i = Ae,(0 + tt), (10) 

where is the average of 6^(0) over 0. An important observation in Fig. 3.5 is that |Aej(0)|'s 
have significant magnitude only on edges (curves Fi and F2) and at corners (curve C). For 
example, |Aej(0)|'s at sites next to an edge (curves Si and S2) are smaller than those of 
curves Fi and F2 by an order of magnitude. Curves I at sites further inside become quickly 
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Fig. 3.5. The local energy at each site i, €,(0), of the 0-LT state, which is a single domain state with 
the LT order of a fixed 0, of the $ = square with L e =16. Curves F, S, C, N and I denote £j(0) 
at sites on the first (green and blue) and second (light green and light blue) lines from edges, at 
the first (black) and the two next (gray) points along the diagonals, and at further interior sites 
(red). The ©-dependence of the total energy £ (0) divided by four is also shown by setting its 
average over (= —143.42) at that of curve C, ec (= —1.347). 



independent of G. Moreover, since Aej(0)'s on sites Fi and F2 nearly cancel with each other 
for a fixed O, the G-dependence of £(0) is expected to come out dominantly from ej(G)'s 
at the four corners. This is in fact the case as seen in Fig. 3.5 where we also draw the 
curve of AS (G)/4 = [£(0) — £]/4. Here £ is the average of £(0) over 0. We emphasize that 
£/4 = -143.42 whereas |A£(0)|/4 ~ [|Aej(0)j on the corner] ~ 0.24 as shown in Fig. 3.5. The 
0-LT state with minimum £(0) is the micro- vortex state with = 37r/4 whose four corner 
spins are in the vortex configuration, while that with maximum £(0) is also the micro-vortex 
state but with = 7r/4 whose four corner spins are in pairs of the d-d spins. 

As compared with the minimum £(0) (= —573.96 with = 37r/4), the ground state 
numerically simulated has the total energy, £ G = —587.33. Its multi-domain structure is 
shown in Fig. 3.4(a-i). The latter is what we can simply think of from the above observation 
on the energetics of the 0-LT states. All spins on the four edges align in parallel to each edge 
they lie on in the simulated ground state. This yields energy gain relative to £(7r/4) by the 
amount estimated as ~ 4 x (L — 2) x [|Aej(0)| of curve Fi or F2] ^ 18.51. This figure is 138% 
of the difference £(tt/4) — £ G . Note that the ground-state energy gain from the four corner 
spins in the micro-vortex pattern is already taken into account by this energy difference. We 
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(a) 




Fig. 3.6. The local energies of the simulated ground state {ep} of the <f> = square with L e =16 (a). 
Those of the boundary sites are omitted since they are out of the scale of the figure. The local 
energy differences of the excited state from the ground states, {ef — ef } (b). The difference at the 
two corners with the d-d spins (~ 0.53) is cut off. 

conclude therefore that the major gain of the ground-state energy relative to those of the 
0-LT states is of the order of L, and can be attributed to the reduction of the symmetry of 
local energies in Eq. (3) from 0(2) in the bulk to Z2 on the edges and corners. 

To go one step further into the energetics of the ground state, let us examine fine details 
of its local energies, {ep}- For example, the value ep of the corner spin is —1.729 as compared 
with —1.588 which is the minimum value of curve C in Fig. 3.5, while the averaged value 
of ep's of the edge spins is —2.247 as compared with —2.252 which is the minimum value of 
curve Fi or F2. These fine details, however, give rise to only several percents modification on 
the above estimation of the ground-state energy gain. The major part of the modification, 
or the excess energy of the ground state relatively to that of the 0-LT state comes out from 
ep's of the interior spins which are shown in Fig. 3.6(a). Note that, in the figure, the scale 
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of its abscissa is enlarged and ep's of the edge and corner spins, which are out of the scale 
of the abscissa, are omitted. In each domains with either x- or y-directed af-FMC order, 
ep's are nearly constant and their values are close to that of curve I in Fig. 3.5. On their 
boundaries, i.e., on the region of domain walls, ep's are larger than those of spins within the 
domains. However, their differences are small and even their sum is at most about one fourth 
of the energy gains due to the boundary spins discussed above. The circumstances are the 
same for spins near the center of the square where the four domain walls meet and the spin 
configuration is close to that of the micro-vortex state (Fig. 3.4(a-i)). Although we cannot 
define accurately the region of the domain walls, we can say that its width is significantly 
smaller than what we have in mind for domain walls in ordinary Heisenberg ferromagnets. 

We are faced to a similarly peculiar aspect of domain walls, when we examine the energetics 
of the excited state shown in Fig. 3.4(a-ii), where the af-FMC order of one of the domains in 
the ground state is reversed. In Fig. 3.6(b), we show differences {ep — ep}, where ef is the 
local energy at site i in the excited state. The corresponding total excited energy is only 1.91, 
for which 1.06 arises from the two d-d spins at both ends of the reversed domain and 0.63 
from the 18 sites neighboring to the corners. The pattern of {ef — ep} at the remaining sites 
is rather flat. This implies that the domain- walls both in the ground state and the excited 
states associate excess energies of an almost equal magnitude, and is naturally expected from 
the periodicity of ej(O) represented by Eq. (10). An important result here is that the excess 
energy of the excited state is almost confined at the two corners with d-d spins and at sites 
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Fig. 3.7. Freezing spin configurations in one of the cooling runs of the <t> = L e = 16 dipolar square. 
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neighboring to them. 

The spin configuration on the L e =15 (odd L e ) dipolar square is intrinsically incompatible 
with the 2x2 periodicity of the LT order. Without going into details of the energetics of its 
low-lying states, we here only point out the following aspects. If we remove spins, say, on the 
top and right edges from the ground-state configuration of the L e =16 square (Fig. 3.4(a-i)), 
we obtain the ground state of the L e =15 square shown in Fig. 3.4(b-i). The same procedure 
on the (a-ii) configuration yields the lowest excited state of (b-ii). The boundary spins draw 
a closed loop in the next excited states but with an apparent point defect at the center as 
seen in (b-iii) and with an apparent line defect next to the right edge in (b-iv). Still, the total 
energies of the two states almost coincide with each other, and are higher than the ground- 
state energy of (b-i) by only 2.24, i.e., of the order of L°. This is another observation on the 
peculiar nature of domain walls as well as defects and dislocations in a finite dipolar system. 

Now let us discuss the freezing characteristics of the dipolar squares. We can see in Fig. 3.3 
that the freezing branches are fixed, or thermally blocked, at around T ~ 0.4 in the present 
simulation. We call this temperature as the domain-blocking temperature and denote it as 
Tqb- In Fig. 3.7, we show the freezing patterns {Si} observed in one of the cooling run of the 
L e = 16 square which reaches to the ground state. These figures lead us to the from-edge-to- 
interior freezing scenario 6 also for the present dipolar square. Once the short-range LT order 
reaches the system size at around T* , at which specific heat exhibits maximum (Fig. 3.1), 
spins on the edges and their vicinity tend to align in an af-FMC order due to the symmetry 
reduction to Z2 on the edges. As the temperature decreases further, each af-FMC domain 
grows, thereby it looks for a proper alignment with other domains while it is overturning 
due to thermal fluctuation. There also appear and disappear domains with the micro-vortex 
LT order (Fig. 1.1(c)) inside the system. The consequence is an almost common increase of 
the freezing parameter S and the order parameter Mlt as seen in Fig. 3.2 and Fig. 3.3 (a), 
respectively. At around T ~ Tdb domains (or spins on the edges) are thermally blocked. At 
lower temperatures than Tdb> S and Mlt further increase simply due to the reduction of 
thermal fluctuations on individual spins within the domains. But Mlt does not saturate to 
unity even at the lowest temperature simply because of the multi-domain structure of the 
ground state associated with the peculiar energetics discussed above. 
3.3 & = 7r/4 dipolar squares 

From the behavior of freezing parameter S shown in Fig. 3.2, one may expect that the 
freezing characteristic of the = 7r/4 dipolar square with = 22 (even L^) having a pair 
of spins at each corner and that of the Ld = 23 (odd L^) square having a single spin at 
each corner are not much different from each other. Actually this is also the case for Mlt's 
observed in eight independent cooling processes, the result of which is shown in Fig. 3.8 only 
for the even-L^ square. In contrast to the $ = square, Mlt's of the <I> = 7r/4 square both 
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Fig. 3.8. The order parameter Mlt of $ = 7r/4 squares with = 22 in each cooling run. 

with even- and odd-L^ simulated end up with one branch at lowest temperatures and their 
values reach to unity within our numerical accuracy. When, however, the spin configurations 
at the lowest temperature are looked at carefully, as shown in Fig. 3.9, essentially only one 
pattern (a) is found in the even-L^ square, while two are found for the odd-L^ square (b-i, 
ii). The three patterns consist of a single domain of the af-FMC order in parallel to one of 
the diagonals of the square (or one of the lattice axes) but with different modulations starting 
from both ends of the other diagonal. 

Let us first introduce the 0-LT states and examine their energetics as in the previous 
subsection. Figure 3.10(a) shows ej(0)'s of the even-L^ square. One can see clearly the short- 
ranged characteristics of Ae;(0)'s also in the present square, namely, Aej(0)'s at interior 
sites only a few lattice distances from the boundary become almost independent of 0. One 
of the characteristic aspects of the figure is that Aej(O) at edge sites (curves F) and at sites 
next to the edges (curves S) have almost equal amplitudes and are almost perfectly out-of 
phase to each others. This implies that their sum (adding Aej(O) at one of the two N sites 
to the sum), which we may regard as the sum of Aej(0)'s at sites on a zigzag edge, becomes 
independent of O. The consequence is the absence of energy gain of the order of L from each 
zigzag edge, in a sharp contrast to the geometry effect in the $ = square. Another aspect to 
be noted is that there appear two types of Cj(O) at the corner sites, curves Ci and C2, whose 
Aei(0)'s have also equal amplitudes and are perfectly out-of phase to each others. Reflecting 
these aspects due to the geometry of the present dipolar square, the total energy of the O-LT 
state turns out to be almost independent of O, i.e., £ = —569.54 with |A£(©)| < 0.86. The 
above arguments can be straightforwardly extended to the energetics of the odd-L^ square, 
by taking into account the further fine detail that both curves F and S are split into two types 
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Fig. 3.9. Spin configurations at the lowest temperature obtained in the cooling processes shown in 
Fig. 3.8: (a) = 22 and (b) = 23. The solid line in (b) indicates a dislocation line. 

as shown in Fig. 3.10(b). We obtain 8 = —569.33 with |A£"(0)| ~ 0.00 for this dipolar square. 

The total energies of the simulated ground state are lower than 8 above described by 2.07 
and 2.17 for the even- and odd-L^ squares, respectively. As seen in Fig. 3.9(a), the ground 
state of the even-L^ square consists of a single domain with the af-FMC order directed in 
parallel to one of the diagonals of the square with dislocations of the O = ir/4 type which 
start from both ends of the other diagonal and go deeply into the interior almost up to the 
center of the square where the two dislocations meet together smoothly. At and near these 
two ends where the micro-vortex units appear, the ground state gains energy, while it looses 
energy at and near the dislocations inside the square relatively to £j(0 = 0). The situation is 
similar but a little complicated for the odd-L^ square. In its ground state (Fig. 3.9(b-i)), there 
appear two d-d spins at two corners of the other diagonal. The dislocations of the = 7r/4 
type start from sites next to the corners and go deeply into the interior almost up to the 
center of the square where the two dislocations meet together but are separated by one lattice 
distance in this case. Our cooling simulations on this square end up also with the excited state 
shown in Fig. 3.9(b-ii). In the excited state, the two d-d spins are antiparallel to each other. 
The centers of the two dislocations lie on an identical lattice axis, and the two dislocations 
are separated significantly with each others by the region with the complete af-FMC order. 
We note that the excess energy of the excited state from the ground-state energy is only 
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Fig. 3.10. The local energies e 4 (9) of the 6-LT states of the $ = tt/4, L d =22 (a) and L d =23 (b) 
squares. Curves F, S, C, N and I denote the energies at sites on the first (green and brown) and 
second lines (light green and fresh) from edges, at the corners (black and blue), at sites next to 
the corners (gray and light blue), and at the interior sites (red). Note that in the odd-L^ square, 
curves Ni and N2 are in phase with curves Ci and C2, respectively. 



0.30. Although we can hardly explain these energy gains and losses with enough quantitative 
accuracy, we consider that they represent the nature of the stiffness energy in the spin order 
due to the dipole-dipole interaction, 20 and will be discussed elsewhere. 

We next argue the freezing property of the even-L^ square demonstrated in Fig. 3.11. At 
T = 0.6, which is a little higher than T* , we can see many short-ranged domains of various 
LT orders, sizes of which are a few tenths of L (a). As the temperature decreases, the sizes 
of the domains grow (b), and at such low temperatures as T ~ 0.3, almost a single domain 
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Fig. 3.11. The directions of the spins on the x-y plane and the local freezing parameter Si (a) at 
T=0.6, (b) 0.4, (c) 0.3, (d) 0.0 for $ = tt/4 (L d =22). 



state of the af-FMC order (c-i) as well as that of the micro- vortex state (c-ii) are observed. 
But the pattern is still fluctuating among such nearly single-domain configurations (mostly 
triggered by the discrete, though small, temperature decreases). This corresponds to the two 
branches of Mlt at temperatures from 0.1 to 0.3 in Fig. 3.8. The upper and lower branches 
correspond to the states shown in Figs. 3.11(c-i) and (c-ii), respectively. We specify the domain 
blocking temperature of this dipolar square as Tdb — 0.1 where the lower branch disappears. 
Below TdB) the present MD simulation ends up always with the ground state pattern of 
(d), or Fig. 3.9(a). We interpret this result as follows. Between the nearly = state of 
Fig. 3.11(c-i) and O = 3n/4 state of (c-ii), there exists a free energy barrier whose height 
is very limited, or of a comparable magnitude with the energy difference of the two states 
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itself. The latter is expected to be also small ( < 1.0) by judging from the energy difference 
between the corresponding 0-LT states. The result also implies the absence of other low-lying 
excited states involving such extended defects as dislocations in this system. The freezing 
process of the odd-L^ square (not shown) is almost the same with that of the even-L^ square 
at higher temperatures than Tdb ( ^ 0.1). In contrast to the even-L^ square, however, our 
cooling simulations on the odd-L^ square end up also with the excited state of Fig. 3.9(b-ii) 
at the lowest temperature, implying that the excited state is separated from the ground state 
of (b-i) by a free energy barrier of a height significantly higher than their energy difference, 
£ — £ G ~ 2.1. These minute but very peculiar aspects of the free energy structure are, we 
consider, the consequence of the intrinsic frustration effect of the dipole-dipole interaction 
combined with the geometry effect of these dipolar squares. The proper understanding of its 
details is remained as a future problem. 

4. Summary and Discussion 

By the MD simulation based on the LLG equation at finite temperatures, we have inves- 
tigated the magnetic ordering and freezing properties of the simplest finite dipolar systems, 
i.e., the dipolar squares cut out from a square lattice. Distinctly different phenomena are 
found depending on how a square is cut out (geometry effect). The ground state of the $ = 
square, whose edges are parallel to the lattice axes, is found to have a multi-domain structure 
which consists of domains having the two types of the bulk LT (af-FMC) order. The peculiar 
"from-edge-to-interior freezing" to the ground state is observed as in the dipolar cube. 6 These 
results are attributed to the reduction of the symmetry of the dipole-dipole interaction energy 
from 0(2) in the bulk to Z2 on the edges. The symmetry reduction is, in turn, attributed to 
the peculiar anisotropic nature of the dipole-dipole interaction that it depends on the rela- 
tive direction of sites in a real space of the interacting spins. In the $ = 7r/4 square, whose 
edges are rotated by 7r/4 from those of the $ = square, the effect of the same symmetry 
reduction as above is canceled in the local energies of each of the four zigzag edges of the 
square. Its ground state consists of a single domain with the af-FMC order in parallel to the 
one of the diagonals of the square with modulations starting from both ends of the other 
diagonal. The freezing starts from the interior, and the final thermal blocking to a state at 
lowest temperatures is governed by spin configurations at and near both ends of the other 
diagonal. 

The excess energy associated with a domain wall in the multi-domain ground state of 
the <I> = square is evaluated to be of the order of L° with L being a linear dimension of 
the square. This is similar to that of a domain wall in an ordinal Heisenberg ferromagnetic 
square under the antiperiodic boundary condition. In contrast to the latter domain wall whose 
width is of the order L, however, that of the present dipolar square seems rather small though 
we cannot specify it accurately. In the ground state of the $ = 7r/4 square, we observe 
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dislocations in the af-FMC order directed in parallel to one of the diagonals of the square. 
They start from both ends of the other diagonal and go deeply into the interior. Still, their 
associated excess energy is of the order of L° and their width looks rather small, as is the case 
for the domain wall in the $ = square. We consider that these characteristics of apparently 
extended objects, domain walls and dislocations, reflect the intrinsic frustration effect of the 
dipole-dipole interaction. The frustration effect, which yields the peculiar ground state with 
the continuous 0(2) symmetry already in the bulk square lattice, makes it possible for spins 
to adjust to further (geometrical) restriction imposed in a finite system only with a quite 
small excess energy. 

In the pattern of the local magnetic energies of the ground states, {ep}, its major changes 
induced by the presence of the boundary are seen to propagate to the interior by only one or 
two lattice distances from the boundary. Their spatial dependence associated with the domain 
walls and dislocations mentioned above is quite small in the sense that the total excess energy 
of an excited state associated with such a defect is only of the order of L°. It is rather surprising 
that the long-ranged nature of the dipole-dipole interaction seems not to play an essential role 
on determining the ground-state and excited-state configurations. To get a deeper insight on 
roles of the long-ranged nature of the interaction, we have preliminarily carried out a similar 
analysis to that we adopted in our previous study on the dipolar cube, 6 namely, repeating the 
identical simulations but with the dipole-dipole interaction whose range is artificially cut off, 
and comparing the results of the two simulations. The consequence of the $ = square is that 
the long-ranged nature makes easier for the system to reach the ground-state multi-domain 
structure as observed in our previous work. For the 3> = tt/4 square, on the other hand, the 
clear consequence has not been obtained yet because of the small energy differences between 
the states with and without dislocations. The problem remains for a future study. 

Lastly we append an explanatory note on our MD simulation. Most of the freezing char- 
acteristics reported in the present paper is rather robust to the cooling rate of the simulation 
under a fixed value of the temperature decrease AT. For example, a rapid cooling with r 
one order smaller than the present one does not yield a significant change either in the peak 
temperature of the specific heat T* or the thermal blocking temperature Tdb- In this sense 
we consider that our choice of the parameters r and AT is appropriate for extracting the 
characteristic properties of the ground state and the low-lying excited states, as well as the 
freezing characteristics of finite dipolar squares of our interest. From a technological point of 
view on nanomagnetism, the value of Tdb of a given dipolar square (whether it is higher than 
the room temperature or not) is certainly of importance. To answer this question, we need to 
know, not only the systematic L-dependence of the energy difference between the ground and 
excited states that we have investigated in the present work, but also that of the free energy 
barrier height between the states. But such an analysis is beyond the scope of the present 
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